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5^ Abstract. We consider viscosity solutions to nonlinear uniformly parabolic equations in 

^J^ nondivergence form on a Riemannian manifold M, with the sectional curvature bounded 

.^^ from below by -k for a: > 0. In the elliptic case, Wang and Zhang | WZl recently extended 

the results of ICal to nonlinear elliptic equations in nondivergence form on such M, where 
^^ they obtained the Hamack inequality for C^- classical solutions. We establish Harnack 

fSJ inequality for nonnegative viscosity solutions to nonlinear parabolic equations in nondi- 

vergence form on M. Harnack inequality of nonnegative viscosity solutions to the elliptic 
equations is also proved. 
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1 . Introduction and main results 



In this paper, we study Harnack inequality of viscosity solutions to nonlinear uniformly 
parabolic equations in nondivergence form on Riemannian manifolds. Let {M,g) be a 
smooth, complete Riemannian manifold of dimension n. We consider a nonlinear uni- 
formly parabolic equation 

(1) F{D^u)-d,u^f in MxR, 

where D^u denotes the Hessian of the function u defined by 

for any vector fields X, Y on M, and Vm is the gradient of u. We notice that in the case, 
when F is Laplacian, ([T]i is the well-known heat equation with a source term. 

In the setting of elliptic equations on a Riemannian manifold M, Cabre llCall established 
Krylov-Safonov type Harnack inequality of classical solutions to linear, uniformly elliptic 
equations in nondivergence form, when M has nonnegative sectional curvature. Krylov- 
Safonov Harnack inequality is based on the Aleksandrov-Bakelman-Pucci (ABP) estimate, 
which is proved using afline functions in the EucUdean case. Since affine functions can not 
be generalized into an intrinsic notion on Riemannian manifolds, Cabre considered the 
functions of the squared distance instead of the affine functions to overcome the difficulty. 
Later, Kim 21 improved Cabre's result removing the sectional curvature assumption and 
imposing the certain condition on the distance function (see fK" p. 283]). Recently, Wang 
and Zhang |WZ| obtained a version of ABP estimate on M with Ricci curvature bounded 
from below, and Harnack inequality of classical solutions for nonlinear uniformly elliptic 
operators on M with a lower bound of the sectional curvature. 

In the parabolic case, Kiylov-Safonov Harnack inequality was proved in PKKL| for 
C^'- classical solutions to linear, uniformly parabolic equations in nondivergence form, as- 
suming essentially the same condition introduced by Kim ||Kl. ABP-Krylov-Tso estimate 
discovered by Krylov |Kr| in the Euclidean case (see also EIWl) is a parabolic analogue 
of the ABP estimate, and a key ingredient in proving paraboUc Harnack inequality. In 
order to prove ABP-Krylov-Tso type estimate on Riemannian manifolds , an intrinsically 
geometric version of Krylov-Tso normal map, namely, 

0(x,t) :- lexp^'V ^u(x,t), — cfi {x,e\p^'Vu(x,t)) - u(x,t)\ 



was introduced. The map <J> is called the parabolic normal map related to u{x, t) and the 
Jacobian determinant of <1> was explicitly computed in |KKL Lemma 3.1]. 

In this paper, we shall prove Krylov-Safonov Harnack inequality of a class of viscosity 
solutions to ([T]l on M with the sectional curvature bounded from below. Let Sym TM be 
the bundle of symmetric 2-tensors over M. A nonlinear operator F : Sym TM — » R will be 
always assumed in this article to satisfy the following basic hypotheses: 

(HI) F is uniformly elliptic with the so-called ellipticity constants Q < A < A, i.e., for 
any S € Sym TM, and for any positive semidefinite P e Sym TM, 

A trace(P.,) < F(S ,, + P.,) - F(S ,,) < A trace(P.J, Vx e M. 
We also assume that F(0) = 0. 
We may assume that < A < I < A. Now we recall viscosity solutions, which are proper 
weak solutions for nonlinear equations in nondivergence form. In the Euclidean case, 
existence, uniqueness and regularity theory for the viscosity solutions have been developed 
by many authors (see for instance, IICILIICCllWll ). In BAFSI IZll. the concept of viscosity 
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solutions has been naturally extended on Riemannian manifolds, which can be found in 



Definitions 2.13 and 5.3 The authors in lAFSIlPZllZl have shown comparison, uniqueness 
and existence results for the viscosity solutions on Riemannian manifolds. 

To obtain Harnack inequality of viscosity solutions from a priori estimates in Subsec- 
tion 4.1 IKKLI . and OCal iKl IWZI . we use regularization by sup and inf-convolutions, 
introduced by Jensen |T|. The classical ABP estimate for viscosity solutions was proved 
by making use of affine functions, especially the convex envelope of the viscosity solu- 
tion (see ICCIIWl ). Replacing affine functions by the squared distance functions on M, as 
mentioned above, we consider the sup and inf- convolutions on Riemanninan manifolds 
defined as follows: For e > 0, the inf- convolution of u, denoted by «£, is defined as 

Us(x) := inf \ u(y) H d^(y, x)} , Vjc e Q c M. 

yen I 2s J 

The sup- convolution can be defined in a similar way using concave paraboloids. We 
show that the regularized functions by the sup and inf- convolutions are semi-convex and 
semi-concave, respectively, which will imply that they admit the Hessian almost every- 



where thanks to Aleksandrov theorem (|Al|Bl. It is proved in Lemma 3.6 that regularized 
viscosity solutions solve approximated equations in the viscosity sense, provided that the 
sectional curvature of M has a lower bound and the operator F is intrinsically uniformly 
continuous with respect to x; see Definition |3.4| As in the Euclidean case, we shall con- 
sider a class of all uniformly elliptic operators, which can be expressed in terms of Pucci's 



extremal operators; see Definition 2. 14 and the comment after it. So, intrinsic uniform con- 
tinuity of Pucci's operators suffices to obtain uniform estimates for the uniformly elliptic 
operators since a class of all viscosity solutions for uniformly elliptic operators is invariant 
under the regularization processes of sup and inf- convolutions. Thus the application of a 
priori estimates to the sup and inf- convolutions of viscosity solutions will give Harnack 
inequality of viscosity solutions. 

On the other hand, assuming the sectional curvature of M to be bounded from below. 



we establish the parabolic Harnack inequality of classical solutions in Section 4. 1 influ 



enced by Wang and Zhang | WZ|, who studied the elliptic case. We introduce the parabolic 



contact set ^a,h for a,b > in Definition 4.3 which consists of a point (x, t) e M xR, 
where a concave paraboloid 

a 9 
--d^Xx) + bt + C (for some C) 

touches u from below at (x, f) in a parabolic neighborhood of {x,l), i.e, Br(x) x (f - r^,~t] 
for some r > 0. Under the assumption that the sectional curvature of M is bounded from 
below, an estimation of the Jacobian of the parabolic normal map on the parabolic contact 



set ^a,b is essential to prove a prior Harnack estimate. Proposition 4.8 



Now we state our main results as follows. In the statements and hereafter, we denote 



Jq^' \q\Jq^ 



and 



K,-(xo, to) := Br(xo) x (to - r^, to], (xq, to) e M x\ 



where \Q\ stands for the volume of a set 2 of M or M x R, and Br(xo) is a geodesic ball of 
radius r centered at xq. 

Theorem 1.1 (Parabolic Harnack inequality). Assume that M has sectional curvature 
bounded from below by -k for k > 0, i.e.. Sec > —k on M, and F satisfies (HI). Let 
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/ e C (K2RixQ,4R^)) ■ If u e C(K2RixQ,4-R^)) is a nonnegative viscosity solution of the 
equation F{D^u) — d,u — f in K2r{xo,4R^), then we have 



u + R^ij- 



(2) sup u<C\ inf u + R^f 1/1"*'+' 

where 6 :— 1+ logj cosh(8 yficR) and C > is a uniform constant depending only on n. A, A 
and yfkR. 

Theorem 1.2 (Weak Harnack inequality). Assume that Sec > -k on M for k >Q, and F 
satisfies (HI). Let f ^ C (K2Rix(),4R^)] .IfueC (A'2r(xo, 4R^)] is a nonnegative viscosity 
supersolution of the equation F{D^u) — dfU — f in K2Rix^),4R^), then we have 

l-f u"]" <C\ inf ^u+R^ I -f irr'Y^' I ; r := max(/, 0), 

\JKg(xo,2R') I l^KsixoAR-) \Jk2r(xoAR-) I j 

where :— I + log2 cosh(8 ^fKR), and the positive constants p e (0, 1) and C are uniform 
depending only on n. A, A, and yficR. 

In the elliptic setting, we have Harnack inequalities of viscosity solutions. 

Theorem 1.3 (Elliptic Harnack inequality). Assume that Sec > -k on M for k>Q, and F 

satisfies (HI). Let f & C (B2r(xo)) .IfueC (B2r(xo)) is a nonnegative viscosity solution of 
the equation F{D^u) — f in B2r{xq), then we have 

(3) sup M<c| inf M + 7?2(-f |/r)H, 

B«(.vo) \Bs_{xo) XJBjgixa) j j 

where :— \ + log2COsh(8 ^[KR\ and C > Q is a uniform constant depending only on 
n,A,A, and yfkR. 

When A- = 0, (|2| and (|3]l become global Harnack inequalities which extend the classical 
Euclidean theory of Krylov and Safonov |KS|. 

Theorem 1.4 (Weak Harnack inequality). Assume that Sec > -k on M for k >Q, and F 

satisfies (HI). Let f e C (B2r(x(j)) . Ifu e C (-B2r(-^())) is o nonnegative viscosity supersolu- 
tion of the equation F(D^u) = / in B2r(xq), then we have 



Jb„ 



uP\ <C\mi u + rH\ |.rr ; r:=max(/,0), 



B«(-i-o) / (^Sfi(.vo) \Jb2r(ao) 

where Q :— \ + logj cosh(8 V^^) and the positive constants p e (0, 1) and C are uniform 
depending only on n. A, A, and yfi<R. 

The rest of the paper is organized as follows. In Section l2] we give some results on 
Riemannian geometry and viscosity solutions that are used in the paper. In Section l3] we 
investigate basic properties of the sup and inf- convolutions, and the interaction between 
the viscosity solution and its sup and inf- convolutions. Section Hlis devoted to prove 
the parabolic Harnack inequalities of viscosity solutions. In Sectiorip] we prove Harnack 
inequalities of viscosity solutions to the elliptic equations. 
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2. Preliminaries 

2.1. Riemannian geometry. Let {M,g) be a smooth, complete Riemannian manifold of 
dimension n, where g is the Riemannian metric and Vol := Vol^ is the Riemannian measure 
on M. We denote {X, Y) := g{X, Y) and |Xp := {X,X) for X,Y e T^M, where T,M is the 
tangent space at x e M. Let d(-, ■) be the distance function on M. For a given point y e M, 
dy(x) denotes the distance function to y, i.e., dy(x) := d(x,y). 

We recall the exponential map exp : TM — > M. If y , x : M. — > M is the geodesic starting 
at X G M with velocity X e T^M, then the exponential map is defined by 

exp,.(^) := r.,z(l). 

We observe that the geodesic y^i x is defined for all time since M is complete. For X e T^M 
with |X| = 1, we define the cut time tdX) as 

tciX) :- sup jf > : exp^{sX) is minimizing between x and exp^(tX)] . 

The cut locus of x e M, denoted by Cut(x), is defined by 

Cut(x) := {exp^(t,{X)X) : X e T^M with \X\ = 1, t,(X) < +00) . 

If we define 

E^ := {tX e r.iM : < f < tdX), X e T^M with \X\ = 1) c r^M, 

it can be proved that Cut(x) - exp^(dE^),M - exp^iE^) U Cut(x), and exp^ : E^ — > 
exp^(£^) is a diflFeomorphism. We note that Cut(x) is closed and has measure zero. Given 
two points X and y i Cut(x), there exists a unique minimizing geodesic txpJ^tX) (for 
X G E^) joining x to y with y - exp^_(X), and we will write X - exp;^'(y). For any 
X i Cut(y) U [y], the distance function c/,, is smooth at x, and the Gauss lemma implies that 

„, , , exp;^'(y) 

Vdyix) = -■ —— 

|exp/(y)| 

and 

V(4/2)(x) = -exp;'(y). 
The injectivity radius at x of M is defined as 

iM{x) :- sup{r > : exp^ is a diffeomorphism from Br(0) onto Br{x)]. 

We note that iM{x) > for any x G M and the map x i-> iMix) is continuous. 
We recall the Hessian of a C^- function u on M defined as 

D^M(X,y):=<VxVM,y), 

for any vector fields X, Y on M, where V denotes the Riemannian connection of M, and 
Vm is the gradient of u. The Hessian D^u is a symmetric 2-tensor in Sym TM, whose value 
at X G M depends only on u and the values X, Y at x. By a canonical identification of the 
space of symmetric bilinear forms on T^M with the space of symmetric endomorphisms 
of TjcM, the Hessian of m at x G M can be also viewed as a symmetric endomorphism of 
T,M: 

D^u(x) ■ X = VxVm, VX g T^M. 

We will write D^u(x) {X, X) = [o^uix) ■ X, x) for X G T,,M. 

Let f be a vector field along a differentiable curve y : [0, a] — > M. We denote by 
-gj(f) = Vy(,)^(f), the covariant derivative of ^ along y. A vector field ^ along y is said to be 
parallel along y when 

-^it) = on [0,fl]. 
at 
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If y : [0, 1] — > M is a unique minimizing geodesic joining x to y, then for any ^ e T-^M, 
there exists a unique parallel vector field, denoted by L^y[,{t), along y such that L:^y(iQ) - 
f . The parallel transport of ( from xioy , denoted by Lx^y(, is defined as 



which will induce a linear isometry L,^. : Tj^M — > T^.M. We note that L^. ^ - L^]. and 



L,,,^ := L,,,^(l) € r,M, 

,, : r,M ^ r,,M. W. ..... ..... ^3,, - ^,3, 

(4) (L,,/,v)^=(^,L,,,y)^, ^C^TM, v e T,M. 

We also define the parallel transport of a symmetric bilinear form along the unique mini- 
mizing geodesic; see liAFSi p. 311]. 

Definition 2.1. Let x,y e M, and let y : [0, 1] — > M be a unique minimizing geodesic 
joining x to y. For S e Sym TM^, the parallel transport of S from x to y, denoted by 
Lj:,v o S ,is a symmetric bilinear form on T^M satisfying 

((l,,, oS)-v,v)^:^{s ■ (l,, ,v) , L,,,v)^ , Vv e TyM. 

Identifying the space of symmetric bilinear forms on TyM with the space of symmetric 
endomorphisms of TyM, L^-y o S can be considered as a symmetric endomorphism of TyM 
such that 

(l,,v o 5 ) ■ V = L^,,. (5 ■ (iv.A v)) , Vv e TyM. 

Then it is not difficult to check that S and L^y o S have the same eigenvalues. 
Let the Riemannain curvature tensor be defined by 

R(X, Y)Z = VxVyZ - VyVxZ - V[X,F]Z- 

For two linearly independent vectors X,Y e T^M, we define the sectional curvature of the 
plane determined by X and Y as 

{R(X, Y)X, Y) 

Let Ric denote the Ricci curvature tensor defined as follows: for a unit vector X e TyM 
and an orthonormal basis {X,e2,- ■ ■ , e„] of TyM, 

n 

Ric(X,X) = YjSec{X,ej). 

As usual, Ric > Kon M (/c e R) stands for RiCj^ > Kgy for all x e M. 

We recall the first and second variations of the energy function (see for instance, iDll. 

Lemma 2.2 (First and second variations of energy). Let y : [0,1] ^ M be a minimizing 
geodesic, and ^ be a vector field along y. For small e > 0, let h : (-e, e) X [0, 1] -^ M be a 
variation ofy defined as 

h(r,t) := exp^(,) r^(0. 
Define the energy function of the variation 

^' dh 2 



E{r) 



Jo 



8t^'^'^ 



dt, for r e (-e, e). 



Then, we have 
(a) 

E{Q) ^ dHjiQ),y(l)) , 
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(b) 



(c) 



^£'(0) = <^(i),r(i))-<^(0),r(0)>, 

In particular, if a vector field ^ is parallel along y, then we have -§ = Q and (^, y) = C 
(for C e R ) on [0, 1]. In this case, we have the following estimate: 

(5) E{r) = £(0) -'■^ ^ {R (r(f), m) jit), m) dt + o (r^) . 

Now, we state some known results on Riemannian manifolds with a lower bound of the 
curvature. First, we have the following volume doubling property assuming Ricci curvature 
to be bounded from below (see fVl for instance). 

Theorem 2.3 (Bishop-Gromov). Assume that Ric > -(n - \)k on M for /c > 0. For any 
Q < r <R,we have 

/-;^ V0l(B2r(z)) I j /- \ 



We observe that the doubling property (|6| implies that for any < r < R < Rq, 

Vol(B,(z)) - \rl 



where D :- 2" cosh" ' (2^/l<R()) is the so-called doubling constant. Using the volume 
doubling property, it is easy to prove the following lemma. 

Lemma 2.4. Assume that for any z & M and < r < 2Rq, there exists a doubling constant 
£) > such that 

No\{B2r{z))<DNo\{Br{z)). 
Then we have that for any Briy) C Br{z) with < r < R < Rq, 



(7) 



HM^^iLr^r^ -^-- 



In particular, if the sectional curvature of M is bounded from below by —k (k > 0), then (|7]l 
holds with :- I + log2 cosh(4 -s/kRo). 



In the parabolic setting, it follows from Lemma 2.4 that for any < r < R < Rq, 

(8) f \r\fr\ <2a-^A^ \R\f\"'' 

where K,„'-(y, s) := B^iy) x{s- ar^, s] c Kr{z, t) = Br{z) x (f - R^, t] for a > 0. 

We recall semi -concavity of functions on Riemannian manifolds which is a natural gen- 
eralization of concavity. The work of Bangert |B] concerning semi-concave functions 
enables us to deal with functions that are not twice differentiable in the usual sense. 

Definition 2.5. Let Q be an open set ofM. A function (f> : Q. ^ M. is said to be semi-concave 
at xq e Q, if there exist a geodesically convex ball Brixo) with < r < iMixo), and a smooth 
function ^F ; Bi-ixo) — > R- such that cf) + ^' is geodesically concave on Br(xo). A function (p 
is semi-concave on Q. if it is semi-concave at each point in Q. 
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The following local characterization of semi-concavity is quoted from MCMSI Lemma 

3.11]. 

Lemma 2.6. Let (p : Q. ^ 'R be a continuous function and let xq e Q, where D. <z M is 
open. Assume that there exist a neighborhood U o/xq, and a constant C > such that for 
any x € U and X e Tj^M with \X\ - I, 

(f> (exp^. rX) + (f> (exp^. -rX) - 2(p(x) 
hm sup T— < C. 

Then (p is semi-concave at xq. 

Hessian bound for the squared distance function is the following lemma which is proved 
in BCMSI Lemma 3.12] using the formula for the second variation of en ergy. According to 

dl is semi-concave 



2.7 



the local characterization of semi-concavity combined with Lemma 

on a bounded open set Q c M for any y e M, provided that the sectional curvature of M is 

bounded from below. 

Lemma 2.7. Let x,y & M. If Sec > —k (k > 0) along a minimizing geodesic joining x to y, 
then for any X € T,M with \X\ = 1, 

d^.iexp^rX) + d^.iexp,-rX)-2d^Xx) ^ / r \ 

lim sup — '■ — z '■ < 2 yKdy(x) coth I yKdy{x)) . 

The following result from Bangert is an extension of Aleksandrov's second differen- 
tiability theorem that a convex function has second derivatives almost everywhere in the 
EucUdean space E (see also fV) Chapter 14]) . 

Theorem 2.8 (Aleksandrov-Bangert, |B |). Let D. c M be an open set and Zef : Q — > R 
be semi-concave. Then for almost every x e Q, is dijferentiable at x, and there exists a 
symmetric operator A(x) : TxM -^ T^M characterized by any one of the two equivalent 
properties: 

(a) for^e T,M, A(x) ■ £, = V^V0(ji:), 

{b) ,^(exp, ^) = ^{x) + < V0(x), ^) + i {A(x) ■ ^, ^) + o (l^p) as^^ 0. 

The operator A(x) and its associated symmetric bilinear from on T^M are denoted by 
D^ipix) and called the Hessian of(p at x when no confusion is possible. 

Let M and A^ be Riemannian manifolds of dimension n and : M — > A^ be smooth. The 
Jacobian of (p is the absolute value of determinant of the differential d(p, i.e., 

Jac (p(x) :- I det d(p(x)\ for x e M. 

The following is the area formula, which follows easily from the area formula in Euclidean 
space and a partition of unity. 

Lemma 2.9 (Area formula). For any smooth function (^:MxR— >MxR and any 
measurable set E <z M x'R, we have 

I Jac <p(x, t)dV(x, 0=1 'H°[E n 4>-\y, s)]dV(y, s), 

Je JmxR 

where 'H^ is the counting measure. 
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2.2. Viscosity solutions. In this subsection, we consider a refined definition of viscosity 
solutions to paraboUc equations slightly different from the usual definition in IZj; see IIWI 
for the Euclidean case. 

Definition 2.10. Let Q c M be open and T > Q. Let u : Q X (0, T] -^ R be a lower 
semi-continuous function. We say that u has a local minimum at {xq, fo) G f^ X (0, T] in the 
parabolic sense if there exists r > Q such that 

u(x, t) > m(xo, fo) for all {x, t) e K^ixo, to) :— Br{x{)) x (to — r^, fg]. 

Similarly, we can define a local maximum in the parabolic sense. 

Definition 2.11 (Viscosity sub and super- differentials). Let Q. c M be open and T > Q. 
Let u : Q.X (0,T] -^ R be a lower semi-continuous function. We define the second order 
parabolic subjet ofu at (x, t) e D.X (0, T] by 

r^''u{x, t) := [[dMx, t), Vifiix, t), D^ipix, f)) e R x T_,M x Sym TM, : (f e C^'^ (Q x (0, T]} , 
u — tp has a local minimum at (x, t) in the parabolic sense] . 

If {p,(,A) e f^^'^u{x,t), then (p,^) and A are called a first order subdifferential (with 
respect to (t, x)), and a second order subdifferential (with respect to x) ofu at (x, t), respec- 
tively. 

In a similar way, for an upper semi-continuous function m : Q X (0, T] — » R, we define 
the second order parabolic superjet ofu at (x, t) e Q.X (0, T] by 

r^'^uix, t) := [[dMx, t), V(p(x, t), D^ipix, f)) e R x T^M x Sym TMj, : if e C^-' (O x (0, T]) , 
u — (f has a local maximum at (x, t) in the parabolic sense] . 

The following characterization of the parabolic subjet V^-u can be obtained by a simple 
modification of |AFS , Proposition 2.2], |Z, Proposition 2.2]. 

Lemma 2.12. Let m : Q x (0, T] — > R /je a lower semi-continuous function and let (x, t) e 
Q X (0, T]. The following statements are equivalent: 

(a) (p,^,A)e'P'^--u(x,t), 

(b) for^e T,M ando-<0, 

u (exp, f , f + cr) > uix, t) + {(, f ) + cr/7 + - (A ■ ^, ^) + o [\^\^ + |cr|) as (^, a) ^ (0, 0). 

Definition 2.13 (Viscosity solution). LetF : MxRxRxTMxSymTM ^ R,andletD.c 
M be open and T > Q.We say that an upper semi-continuous function m : fi X (0, T] — > R 
is a parabolic viscosity subsolution of the equation dtU = F(x, t, u, Vm, D^u) in Q. x (0, T] if 

p-F(x,t,u(x,t),(,A) <0 

for any (x, t) e Q.X (0, 7"] and (p, (,A) e 'P^'^u(x, t). Similarly, a lower semi-continuous 
function u : D.x(0,T] — > Ris said to be a parabolic viscosity supersolution of the equation 
d,u = F(x, t, u, Vm, D^u) in D. X (0, T] if 

p-F(x,t,u{x,t),^,A) >0 

for any (x, t) € Q.X (0, T] and (p, ^, A) e "P^'^uix, t). We say that u is a parabolic viscosity 
solution ifu is both a parabolic viscosity subsolution and a parabolic viscosity supersolu- 
tion. 
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We remark that parabolic viscosity solutions at the present time will not be influenced 
by what is to happen in the future. In the Euclidean space, Juutinen fJu] showed that a 
refined definition of parabolic viscosity solutions is equivalent to the usual one if compari- 
son principle holds. Whenever we refer to a "viscosity (sub or super) solution" to parabolic 
equations in this paper, we always mean a "parabolic viscosity (sub or super) solution" for 
simplicity. 

We end this subsection by introducing Pucci's extremal operators and their properties. 
We refer to [CC] for the proof. 

Definition 2.14 (Pucci's extremal operator). For Q < A < K (called ellipticity constants), 
the Pucci's extremal operators are defined as follows: for any x e M, and S x & Sym TM^, 

MXa(5.) := M^(Sx) = ^ ^ e,- + A ^ e,, 

e,<() e,>0 

M-,j,(Sx) := M-(Sx) = A ^ e,- + ^ ^ a, 

e,<() i?,>0 

where e,- — e,(5.v-) are the eigenvalues of S x- 

In the special case when A-h- 1, the Pucci's extremal operators M^ simply coincide 
with the trace operator. We notice that the hypothesis (HI) is equivalent to the following: 
for any 5, Pe SymTM, 

(nV) M-(Px)<F{Sx + Px)-F{Sx)<M^{Px\ VxeM, and F(0) = 0. 

Lemma 2.15. Let Sym(«) denote the set ofnxn symmetric matrices. For S,P e Sym(n), 
the followings hold: 

(a) 

M*(S) = sup trace(A5), and M'^(S) = inf trace(A5), 

where Sa,a consists of positive definite symmetric matrices in Sym(n), whose eigen- 
values lie in [A, A]. 

(b) M-(-S)^-M^(S). 

(c) M-(S +P)< M-(S) + M+(P) < M^(S +P)< M+(5) + M^(P). 

Notation. Let r > 0,p > 0, zq & M and to e R. We denote 

K,;p(zo, to) := Br(zo) x (fo - p, h)], 
where B^izo) is a geodesic ball of radius r centered at Zo- In particular, we denote 

Kr{ZO,to) '■- Kr^r^-izoJo)- 

3. Sup and inf- convolutions 

In this section, we study the sup and inf- convolutions introduced by JensenfJl (see also 
BJLSL MCCI Chapter 5]) to regularize continuous viscosity solutions. Let f2 c M be a 
bounded open set, and let m be a continuous function on Q x [To, T2] for T2 > To- For 
e > 0, we define the inf-convolution of u (with respect to Q x (To, T2]), denoted by u^, as 
follows: for (xq, fo) e Q x [T^o, T2], 

(9) UE(xo,to):^ mf \u{y,s) + —{d^(y,xo) + \s - tof-]} . 

(y.s)enx[To,T2] { ^^ J 

Lemma 3.1. For u e C (Q x [To, T2\\ , let u^ be the inf-convolution of u with respect to 
Q X {To, T2\. Let {xo, fo) € Q x [ro, Tj]. 
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(a) IfO < e < s' , then u^'ixo, to) < Ug^xo, to) < u{xo, to). 

(b) There exists (yo, sq) e Qx[ro, T2] such that Ueixo, to) — u{yo, so)+j^ W^(vo,-Xo) + l*o — foR ■ 

(c) cp-(yo, xo) + \so - fol^< 2e\u{xo, to) - u(yo, so)\ < MML''(nx{To,T2])- 

(d) Us 1 u uniformly in D. X [To, T2]. 

(e) Ug is Lipschitz continuous in Q x [To, Ti]: for (jcq, to), {xi, fi) e O x [To, T2\, 

3 3 

(10) \uE{xo,to)-Us(xi,ti)\< — diam(Q)(i(xo,xi)+ — (r2 -ro)|fo-fi|. 

Ze Ze 

Proof. From the definition of u^, (a) and (b) are obvious. From (a) and (b), it follows that 

— [d^iyo, Xo) + \so - fol^j = Usixo, to) - u(yo, so) < u(xo, to) - u(yo, so), 
proving (c). To show (d), we observe that 

< u(xo, to) - Ue(xo, to) < u(xo, to) - u(yo, so)- 

We use (c) and the uniform continuity of m on Q x [To, T2] to deduce that Ug converges to 
u uniformly on O x [To,T2]. 

Now we prove (e). For {y, s) e Q.X [To, T2], we have 

Usixo, to) < u(y, s) + — [d^(y, Xo) + \s - fol^} 

< u(y, s) + Y {("^Cv. -^i) + d(xuxo))^ + {\s - fi | + |fi - fol)^} 

= u(y,s) + — i^d^iy,xi) + d^ixo,xi) +2d{y,xi)d{xo,xi) + (\s - ti\ + |fo-fil)^} 

1 3 3 

< u(y, ^) + ^ [d\y,xi) + \s- tif-] + — diam(Q)£/(xo,xi) + —{T2 - To)\to - fi|. 

Taking the infimum of the right hand side, we conclude ( [TO] i, that is, Us is Lipschitz con- 
tinuous on Q X [To, T2]- □ 

Now, we show the semi-concavity of the inf- convolution, and hence the inf- convolu- 
tion is twice differentiable almost everywhere in the sense of Aleksandrov and Bangert's 
Theoreml2.8l 



Lemma 3.2. Assume that 

Sec > -K on M, for k>0 

For u e C (O x [Tq, T2]] , let u^ be the inf-convolution of u with respect to Q. x (Tq, T2], 
where Q. <Z M is a bounded open set, and To < T2. 

(a) Ue is semi-concave in Q X {To, ^2). Moreover, for almost every (x, t) e D. X (To, T2), 
Us is differentiable at {x, t), and there exists the Hessian D Us(x, t) (in the sense of 
Aleksandrov-Bangert's Theorem \2.8\ such that 

(11) 

Us{s\p^^,t + cr) = Us(x, t) + {Vus(x, t), + crd,Us(x, t) + -{p^Us(x,t) ■^,^\ + o{\^\^ + |o-|) 

as (^, a-) e T,M x R ^ (0, 0). 

(b) D^Us(x,t) < - -v^diam(Q)coth( V^diam(Q)W, a.e. in Q x (To,T2). 

(c) Let H X (Ti, T2] be a subset such that H X [Ti, T2] <z Q.X (To, T2], where H is open, 
and To <T\ < T2. Then, there exist a smooth function ij/ on M X (—00, 7^2] satisfying 

0<\p <\onMx(-oo,T2\, (A s 1 inHx[T^,T2\ and suppiA c Q x (T^o, ^2], 
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and a sequence {w^Y^^i of smooth functions on M X (-00, T2\ such that 

Wic — > i//Ui; uniformly in M X (—00, T2\ as k ^ +00, 

\Vwk\ + \d,Wk\ <C inMx (-00, T2I 

dfWk — » d,Us a.e. in H X (Ti, T2) as k ^ +00, 

D^Wk < Cg in M X (-00, Tt], 

D^Wk — > D^Ue a.e. in H X {T\, T2) as k ^> +00, 

where the constant C > is independent ofk. 

Proof. To prove semi-concavity of m^ in f2 x (To, T2), we fix (xq, to) e Q.X (Tq, T2), and 
find (yo. So) e Q x [To, T2] satisfying 



1 , 2 

ujxo, to) = u(yo. So) + :^\d (jo, xo) + \so - fol 



2e 

For any ^ e Tx„M with |^| = 1, and for small r G R, it follows from the definition of the 
inf- convolution u^ that 

Me (exp^^ r^, to + r) + Me (exp^.^ -r^, to - r) - 2ue(xo, to) 
< u(yo, so) + — [d^ {yo, exp,.^ r^) + \so- (to + '•)l^} 

+ u(yo. So) + — [d^ {yo, exp^.„ -r^) + \so - (to - r)\^} - 2ue(xo, to) 

Then, we use Lemma [ZTJ to obtain that for any ^ e T^^M with |^| = 1, 

Ue {^Wxo '■^' to + r] + Ue (exp^^ -r^, fo - '") - '2.Ue(xo, to) 
lim sup r 

1 < (exp,„ r^) + < (exp,„ -r^) - 2flf2^(^o) 1 

,,,. <limsup— ■ : + - 

(12) ,.^0 2e r^ e 



- V^t/yo(xo)coth(V^(;y„(;i:o)) + - 
< - V'i'diam(Q)coth( y/Kdmm(Q.)) H — , 
where we note that TCOth(T) is nondecreasing with respe ct to t > 0. We recall that Ue is 



< 

s 



Lipschitz continuous on Q x [Tq, T2] according to Lemma 3.1 Since (xq, to) e Qx (To, T2) 
is arbitrary, ( [T2] i and Lemma [Z6] imply that Ue is semi-concave on Q x (Tq, T2). Thus, m^ 
admits the Hessian almost everywhere in Q x (Tq, T2) satisfying ( fTTj ) from Aleksandrov 
and Bangert's Theorem 2.8 The upper bound of the Hessian in (b) follows from ( fTTj l and 

We use a standard mollification and a partition of unity to approximate iJ/Ue by a se- 
quence {wj:)^j of smooth functions in (c), where a mollifier is supported in (-6,0] with 
respect to time (for small S > 0), not in (-6,6). By using Lipschitz continuity of m^ on 
Q. X [Jo, T2] and semi-concavity on Q x (Tq, T2), it is not difficult to prove the properties 
of W/t. For the details, we refer to the proof of Lemma 5.3 in [Ca|. D 

Next, we shall prove that if m is a viscosity supersolution of ([1]), then the inf- convolution 
Me is still a viscosity supersolution; see BCILI Lemma A. 5] for the Euclidean case. 
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Proposition 3.3. Assume that 

Sec > -K on M, for k>0. 

Let H and Q be bounded open sets in M such that H <Z Q., and Tq < Ti < T2. Let 
u e CIQx [To, T2]) , and let oj denote a modulus of continuity of u on Q.x[T^),T2], which is 
nondecreasing on (0, +00) with w(0+) = 0. For e > 0, let u^ be the inf-convolution ofu with 
respect to Q.x(T(), T2\. Then, there exists sq > Q depending only on l|M||/^».(Qxrr t iV^'^' -^"' 
and Ti, such that if < s < eq, then the following statements hold: Let (xQ,to) € H X 
[Ti, T2], and let (yo, sq) e Q.X [Tq, T2] satisfy 

Ue(xo, to) = u(yo, So) + Y {^^Cvo, xo) + ko - fol^} • 

(a) We have that 

(yo,so)enx(To,T2], 
and there is a unique minimizing geodesic joining xq to yo- 

(b) Ifip, (,A) e 'P^'^Ug.ixo, to), then we have 

yo = exp^^(-£(), and sq e [fo - sp, Ti]. 

(c) Ifip, (,A) e P^'^Ugixo, to), then we have 

[P' L.oM^, ivo.vo °A- 2ku\2 ^e||M|li»(nx[7-„,7-,]))^vo) e P^''u(yo, so), 

where Lx„j„ stands for the parallel transport along the unique minimizing geodesic 
joining xo to yo = exp^^(-e^). 



Proof. By recalling Lemma 3.1 (c), we see that 

(yo, so) e B2 ViSi(-^o) x ([fo - 2 ^/^, to + 2 ^/^] n [T^o, T2]) , 

form :- I|m||^»('qx|7- ToiV ^^^^^ '^^e distance between H and 5Q is positive, we select eo > 
so small that 

2 y/mso < min { d(H, dD.), Ti - To\ -: 60, 

where d(H, dQ) means the distance between H and dQ.. For < s < so,we have that 

(yo, so) e fi X (To, T2] 
since (xo, to) e H x[Ti, T2]. We observe that 

i^ :- inf iM(x) : jc € Q } > 

since Q is compact from Hopf- Rinow Theorem and the map x i-> iM(x) is continuous. 
Now, we select 

^o^-^ju^— -minfd?,,/^. 

°ll"llL"(nx[7'o,r2]) 

Then we have that for < e < eo, 

d\xo,yo) < 4e||M|li„(n^[j.„_j.,]) < 4eol|M|li»(nx[r„,r2]) < '|' 

and hence d(xo,yo) < i^ < niin{iM(xo), iM(y())} , which impUes the uniqueness of a mini- 
mizing geodesic joining xo to yo. This finishes the proof of (a). 
From (a), there exists a unique vector X € r^„M such that 

yo = exp^^X, and \X\ = d(xo,yo). 
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First, we claim that if {p,^,A) e P^'^UeixoJo), then yo - exp^.^(-s(), namely, X = -e(. 
Since (p,(,A) e 'P^-'Ui;(xo,t()), we have that for any ^ e T^^M with |^| = 1, small r e 
R, o- < and for any (y, s)eD.x [Tq, T2], 
(13) 

u(y, ^) + ^ {d^ {V' exp,.„ r^) + \s- (to + o-)f} 

> Ue (exp^^ r^, fo + a) 

> M,(xo, fo) + r{^, ^> + £^/5 + ' ^^2^ +o[r^ + \o-\) 

= u{yo, so) + — [cfiyo, xq) + \so- fol^} + r{^, + o-p + ^^^ + o[r^- + |cr|) . 

When (y, 5) = (yo, ^0) and cr = in ( [T3] l, we see that for small r > 0, 
— {^(yo, Jco) + rf > —S (yo, exp^.^ r^) 



2e 2e 

1 



>lrf2(;.o,3'o) + K^,^) + ^^^^^^+o(r^), 



and hence for small r > 0, 

(14) rd{xo,yo)>r{st^,^) + o{P), V^ e r,„M with |^| = 1. 

If X = 0, ([14ji implies that <ef,^) = for all £, e r,„M. Thus we deduce that f = and 
yo = exp^^ = exp,^(-e^). 

Now, we assume that X + ^.\i {y, s) - {yo, so), o" = 0, and ^ = XI\X\ - X/d(xQ,yo) in 
( [T3| ), then we have that for small r > 0, 

^ {fl'(xo,3'o) - rf > ^dHxo,yo) + r{(,0 + !li^_M + ^ (r^) 
and hence for small r > 0, 

(15) -rd(xo,yo)>r{Ei,X/\X\} + 0{r^). 
For small r > 0, ([T4| and ([T5| imply that 

{-e(,^}<\X\^d{xo,yo), "^^ e T,„M with|^| = l, 

and 

{-s(,X/\X\)^\X\=d{xo,yo). 
Then, it follows that -e( = X and hence yo = exp^^ X = exp^_^(-e^) for X i^ 0. Thus we 
have proved that yo = exp^^(-e^). 

When (y, s) - (yo, sq) and r = in ( [T3| l, we have that for small o" < 0, 

1 9 1 2 

— |i0 - k) - cr|- > —\so- to\ +crp + o(|cr|), 
2e Is 

which impUes that so > to - ep. This proves (b). 

To show (c), we recall that there is a unique minimizing geodesic joining xo to yo, 

and (yo, io) £ ^ x (^o, T2] according to (a). Using the parallel transport, we rewrite ( [T3] l as 

follows: for any v e Ty^M with |v| = 1, and small r e R, cr < 0, and for (y, s) e Qx [Fq, ^2], 

^2 
M(y, s) > M(yo, so) + r (^, i>.o,x„ v)^ + crp + — {a ■ {Ly^,,,v) , Ly„,,,v) ^ 



+ y- [d^iyo, xq) - d^ [y, exp,^ rLv„,.,„v)} + — [\so -to\^ -\s-to- o-\^} + o{r^ + |cr|). 
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By setting (y,s) :- (exp,,^ rv, sq + en for small r e R, cr < 0, we claim that 

u (exp,,^ rv, 50 + cr) > u(ya, so) + r (Avo.vo^' ^),,^ + o"/' + y ((Avo.vo ° a) ■ v, v)^^ 

+ :r- {d^iyo, xq) - cf (exp,,^ rv, exp^^ -Ly„,.,„rv)} + o(f + \(r\) 
(16) ^"^ 2 

> uiyo, 5o) + r (l^-„,v„^, v)^^ + crp + — ((l,„,,,„ o a) ■ v, v)^^ 

- —r^Kd^(xo,yQ) + o[r^ + \o-\) . 

The first inequality is immediate from Q and Definition |2.1| To prove the second inequal- 
ity in ( fTS] ), we consider a unique minimizing geodesic 

7(0 := exp^„(-feO 

joining y(0) - xq to y(l) = yo = 6xpj.^(-e^). For a given v e ry„M with |v| = 1, define a 
variational field 

v(0 := Ly^,y(t)y e r^(;)^ 
along 7, where v(0) = i^vo.AoT^; and v(l) - v. For small e > 0, we define a variation /i ; 
(-e,e)x[0, 1] ^M,ofr, ' 

h(r,t) := exp^(,) rv(0. 
The energy is defined as 



Jo 



E{r) 

Jo 






2 



£/f. 



We use the second variation of energy formula (|5]l to obtain 

E{r) = £(0) - '■' J (-^ (r(f), v(f)) r(0, v(0> ^ff + o (r^) 

since y is a unique minimizing geodesic, and v{t) is parallel transported along y. Since 

|v(f)| = |v| = 1, and \y{t)\ = |7(0)| = (^(xcjo) for t e [0, 1], we have that 

£(0) - E(r) ^'■'- [ (R (r(0, v(0) r(f), v(0> dt + o (r^) 

= r^ J Sec (r(0, v(0) ■ (lr(Ol' - <r(0, y(t)f) dt + o (r^) 

> -'•' r ^(lr(Ol' - <r(0, v(0)') dt + o (r^) 

> -r^^ |r(0)|2 + o (r^) = -r^Kd\xo, yo) + o (r^) . 
Recalling that £(0) = d^(xo,y{)), and 

E{r) > d^ (expy(O) rv(0), exp^^j rv(l)) = d^ (exp^.^ A,„,.,„rv, exp^^ rv) , 
we obtain 

d^(xo,yo) - d^ (exp^^ -Ly^^^^rv, exp^^ rv) > £(0) - E(r) 

> -p-Kd^(xQ,y(y) + o{p-j, 
which proves the second inequality of ([16}. 



3.1 



(c), it follows that 
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Since d^{xo,y{)) + \to - sq^ < ^B\\u\\^„r^^,j j,\ from Lemma 

(17) cf-{xQ,yi)) < 2s\u{xi),ti)) - mCvo, io)| < 2s uh. Je|l«ll£»/nx[ro,r,])) ' 

where w is a modulus of continuity of m on Q x [To, Ti]- Therefore, we use ( fTSj l and ( fTTj i 
to conclude that for any v e Ty^M with |v| = 1, and for small r € R, cr < 0, 

u (exp^^ rv, io + cr) > u{ya, sq) + r (Lv„,yo^, v)^^ + crp + y ((Lxo,v„ ° A) ■ v, v)^^ 

„2. 



^/f w (2 ^e||M|li„(nx[7-„.7-,])) + 0(^2 + kl) . 



Therefore, Lemma 2.12 implies 

[p^ Avo,y„^> ^Aojo ° A - 2/f w^2 ^e||M|lz,»(nx[r„,r,])) ^vo j e V^'^uiyo, sq). 



Now, we recall the intrinsic uniform continuity of the operator with respect to x from 
llAFSJ . which is a natural extension of the Euclidean notion of uniform continuity of the 
operator with respect to x. 

Definition 3.4. The operator F : Sym TM — > R is said to be intrinsically uniformly 
continuous with respect to x if there exists a modulus of continuity ojp : [0, +00) -^ [0, +00) 
with a>f{0+) — such that 

(H2) F(S)-F (l,,,, oS)<cof (d(x, y)) 

for any S e Sym TM^, and x,y e M with d(x,y) < min {iuix), iiuiy)] ■ 

We may assume that aif is nondecreasing on (0, +00). Recall some examples of the 
intrinsically uniformly continuous operator from BAFSI . 



Remark 3.5. (a) When M = W,we have L^,,. o S = S so (|H2]l holds. 

(b) In general, we consider the operator F, which depends only on the eigenvalues of 
S e Sym TM, of the form : 

(18) F(5 ) = G ( eigenvalues of 5 ) forsomeG. 

Since S and L^^y o S have the same eigenvalues, the operator F satisfies intrinsic 
uniform continuity with respect to x (with cup = 0). The trace and determinant of S 
are typical examples of the operator satisfying ( |18[ ). 

(c) Pucci's extremal operators AV" satisfy (|18|l, (|H2|l and (HI). 



Lemma 3.6. Under the same assumption as Proposition 3.3 we also assume that F sat- 
isfies (HI) and (|H2}. For / e C (Q x (To, T2]), let u e C(0 x [T^o, ^2]) be a viscosity 
supersolution of 

F(D'^u)-d,u=f inQ.x(To,T2]. 

If < s < So, then the inf-convolution Us (with respect to CI x (To, T2]) is a viscosity 
supersolution of 

F(D\,) - d,u, = /, onHx (r, , Tz], 
where so > is the constant as in Proposition\3.3\ and 



fe(x, t) :- sup f + a)p\2 ■\/mE) + 2nA/c w (2 V^Je) 
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form :- l|M||^»/Qxrr t iV Moreover, we have 

F(D Ui;) - dfUg < fg a.e. in H X {Ti, Ti). 

Proof. Fix < e < eo- Let ip e C^"' {H x (T\ , T2]) be a function such that u^-ip has a local 
minimum at (jcq, Iq) € H x (Ti,T2] in the parabolic sense. Then we have 

[d,(fi(xo, to), V(^(xo, to), D^ipixo, f())) e P^'^Ui^ixo, to). 

We apply Pi'oposition |3.3| to have that 

(d,(p{xo, to), L^„,,,„ V(^(x(), to), L.to.vo ° D^vixo, to) - 2kcl> (l yfmsj §,.„) e P^'^uiyo, sq) 

for 

yo := exp,^ (-eV(^(xo, fo)) e Bj v^(-^o) c Q, 

and some sq e f - 2 V^ne^min jfo + 2 V^ie, 7"2} c (Tq, r2]. Since m is a viscosity super- 
solution in Q X (To, 7^2], we see that 

fiyo, So) > F (l^.„,,.o ° D^(p(xo, to) -2Ka>{2 yjms) g,,„) - d,ip{xo, to) 

^ F (^.vo.vo ° D^fixo, to)) -nK-2K(x)(2 yfms) - d,(f(xo, to) 

> F (D ip(xo, to)] - Ljp id{xo,yo)) - 2nA/<- w (2 y[me\ - d,(p(xo, to) 

using the uniform ellipticity and intrinsic uniform continuity of F. Thus, we deduce that 

F {p^ifixo, fo)) - d,ip{xa, to) < fiyo, so) + (^f id{xo,yo)) + 2nA/^ w (2 Vme) 

< fs(xo,to). 

Therefore, u^ is a viscosity supersolution of F(D^Ug) - d,Ui, = /^ in // x (Ti, T2}. 



According to Lemma 3.2 Ug admits the Hessian almost everywhere in Q x (To,T2) 
satisfying ( [TT| i. For almost every (x, t) e Q x (To, 72), we use ( [TT| l and Lemma 2.12 to 
deduce 

(d,Ue(x, t), WUsix, t), D^Ueix, f)) € P^'^ Us(x, t) H P^'^Ugix, t). 

Therefore, we conclude that 

F{D^Ug) - d,iis < fs a.e. in H x {Ti,T2), 
since m^ is a viscosity supersolution in // x (T^, T2\. □ 



For a viscosity subsolution, we can obtain similar results to Lemmas 3.1|3.2 3.6 and 



Proposition 3.3 using the sup-convolution: 



M''(xo,fo):= sup \u{y,s)- — [cf-{y,xo) + \s-tQ\^Y> fox {xo,to) eQ.x[TQ,T2\, 
0',i)enx[ro,r2] I ^^ > 

under the assumption that the sectional curvature is bounded from below. 

4. Parabolic Harnack inequality 



4.1. A priori estimate. In this subsection, we shall prove Proposition 4.8 which is a main 
ingredient of a priori Harnack estimate. We begin with the definition of the contact set for 
the elliptic case from I. WZJ . 
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Definition 4.1. Let Q.be a bounded open set in M and let u € C(0). For a given a > and 
a compact set E <Z M, the contact set associated with u of opening a with vertex set E is 
defined by 

Ma{E;il;u) := Ix e Q : 3y e E s.t. inf [m + -d^] - u(x) + -d^Xx) 

The following lemma is quoted from BWZI Proof of Theorem 1 .2] and MCal Proof of 
Lemma 4.1] (see also IICMSI Proposition 2.5] and El Chapter 14]). 

Lemma 4.2. Assume that 

Ric > -K on M, for k>0. 

Let Q. be a bounded open set in M and E be a compact set in M. For a > and a smooth 
function u on D., we define the map cf> : Q. ^ M as 

4)(x) '■- exp^a"'VM(x). 

Then, we have the following : Let x e JiaiE; Q; u). 

(a) Ify G E satisfies 

I a r,\ fl 9 

inf M + -dt - u(x) + -drXx), 
n \ 2 ■ / 2 ' 

theny — (j){x) — Qxp^a^^Vu{x), ;c ^ Cut(3'), and -Vu(x) = -dy(x)Vdy(x). 
(b) 

(19) Jac^(.).^"(/^)(^(/^).^" 

\ y n a ) [ \ y n a J na 

where 

^(r) = rcoth(T), ^(t) = sinh(T)/T, t > 0. 

Now we define a parabolic version of the contact set which contains a point (x, 1) 6 

a 9 
M X R, where a concave paraboloid --dy(x) + bt + C (for some a,b > and C) touches u 

from below at (x, 1) in a parabolic neighborhood of (x, 1), i.e., in Kr(x, 1) for some r > Q. 

Definition 4.3. Let Q. be a bounded open set in M and let u € C(Q x (0, T]) for T > Q. 
For given a,b > and a compact set E <Z M, the parabolic contact set associated with u is 
defined by 

y[«,i(£;Qx(0,r];M) 
:= i (x, r) € Q X (0, T]:3yeE s.t. inf (u(z, r) + '^dj(z) - br) = u(x, t) + '^d^ix) -bt\. 

As in KKLl . for u e C^-' {Q. x (0, T]) , we define the map (p : Qx{Q,T] ^ Mhy 

if>(x,t) :- exp^. a"'VM(jc, f), 
and define the parabolic normal map <5 : Q x (0, T] — > M x R by 

(b(x,t) :- icf>(x,t),--d^ (x,(f)(x,t)) - a'^ {u(x,t)-bt} 

Lemma 4.4. Assume that 

Ric > -K on M, for k>0. 
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Let D. be a bounded open set in M, and let u be a smooth function on D. X (0, T] for 
T > 0. For any compact set E <Z M, a,b > 0, and < A < I, we have that if (x, t) e 
J{a.b{E;^x{Q,T\,u),then 
(20) 



Jac <l)(x, t) < 



1 



in + D" 



_^n 



/c|Vm| 
n a 



AAu — dfU b /a^|Vm| 

; + :; — \- n,Jir \ ,/ 

Aa Aa {yna 



nn+l 



where 



^(T) = TCOth(T), J^(t) = sinh(T)/T, T>0. 



Proof Let (x, t) e Jl„_t := ^a,b(E; ^ x (0, T]; u) c Q.X (0, T]. From the definition of the 
paraboHc contact set, there exists a vertex y e E such that 



nx(o,t] 
According to Lemma 4.2 we have that 



inf u{z, t) + -d^.(z) -br]- u(x, t) + -df.(x) ■ 
:x{0,t] \ 2 - / 2 - 



■bt. 



y - <p{x,t) - expj^a Vm(x, f), xiCut(y), and -VuixJ) - -dy(x)'Vdy(x) 



since 



inf I m(z, + -d^Xz)j - u(x, t) + -dy(x). 

We notice that D^ (u + |(i^) (x, t)>0 and d,u(x, t)-b <0. Now we set 

:= </>(■, : Q 3 z H-> exp,.fl"'VM(z, e M 
to obtain from Lemma l4~2l that 

(21) Jac0(x) < ,y^(f-mU([lm^^X\,jy 

\ y n a I [ \ynaj na ] 

By a simple calculation, we have that for (^, cr) e T^M x R\{(0, 0)), 
t/O(A;,f)-(f,o-) = Li0-^ + o--;^, -iv[dll2){y),d^-^ + a--Pi - a'^ o- {d,u - b) , 

where ^(x,0 = ^|^^Q0(x,f + t) e r,,M and we used V{d^.l2){x) = -a-^Vu(x,t). To 
compute the Jacobian of O, we introduce an orthonormal basis {ei , ■ ■ ■ , e„) of Tj^^M and an 
orthonormal basis {ei , ■ ■ ■ , e„] of TyM = T^(x,i)M. By setting for /, j = 1, ■ ■ ■ , n, 

A;y := (e,-, d^ ■ ej) , bi := U, -^\, and c; := (e;, V (^2/2) (y)) , 

the Jacobian matrix of O at (x, t) is 



-Cj^Ajij —c^bii + a (b — d,u) 
Using the row operations and (|2T|l, we deduce that 



Jac (i>(x, t) = 



a (b - o,u) 



— a (b — dtu) Jac ^(x) 



<a-^ib-d,u)y"\f-\^]UI/-\^]^^\(xJ), 
n a \ \ \ n a na 
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where we note that (b - d,u) (x, t) > Q and Jac 0(x) > 0. According to the geometric and 
arithmetic means inequality, we conclude that 



Jac <l)(x, t) < 



1 



(n + If 

1 
(n + 1)" 



1 
(n + 1)" 






k\Vu\ 



^ 



n a 
/<• I Vm| \ \ Xau — dfU + b 



K |Vm|\ AmI b — d,ii 
n a j na \ a 



n+l 



n a 



1 



.n-..r 



Aa 

K |Vm|\1 b - d,u 



n«+l 



J^ 



n a /) a 
K |Vm|\ I Aau — d,u + b 



n a 



Aa 



+ n^e 



a:|Vm| 
n a 



a:|Vm| 
n a 



since (b - d,u) {x, f) > and yij) - sinh(T)/T > 1 > 1 for all r > 0. n 

Assuming the sectional curvature of M to be bounded from below, we have ABP- 
Krylov-Tso type estimate in the following lemma, which will play a key role to estimate 



sublevel sets of u in Proposition 4.8 
Lemma 4.5. Assume that 

Sec > -K on M, for k>0. 

Let Rq > and < rj < I. For zo & M, and < R < Rq, let u be a smooth function in 

KaiR,a^R2(zo,0) c M X R such that 

(22) M>0 in K^^ii^^^R2{zo,Q)\K/j^R^p^R2{zQ,Q) and inf m < 1, 

where ai := — , 02 '■— ^ + jf + \, Pi :— -, and /32 :— 4 + rj^. Then we have 



(23) 

\br{zq)\-r^< r 

J\u<M„\nK,, „ ,, „i 



yn 



'lilK.fyR- 



(20,0) 



III l^"^'^'") - ^'"l +Y^+in + l) j-jA 



where the constant Mjj > depends only on rj > 0, and 
for y{T) = sinh(T)/T, and Jf (r) = TCOth(T). 



(^0,0) 2R 



t\ 



PiR aiR 







1 1 

r,R 


m' 


^2R' 








]rfR' 











aiR' 



Figure 1. a^ := if, az := 4 + //^ + \, p, := }. /?2 := 4 + rf 
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Proof. We consider the parabolic contact set 

^fl,fc(Bfi(zo);^a,«,a2«2feo,0);M) for fl := ^ and b := -^, 

which will be denoted by Ji for simplicity. As in the proof of BKKLI Lemma 3.2], for any 
y G Bgizo), we define 

1 2 1 , b 6 

Wy(x, t) := -R u(x, t) + -dy(x) - C^t, C,, := - = ^. 

From the assumption ( [22] ), we see that 

inf Wv <(5 + -r|/?^ =: A„/^^ 
K2r(zo,0) ■ \ rp-j 



Wy > l6 + ^ j /?2 = (A, + \)R^ on /:„,«, „,R2(zo, 0)\Kp^R,fi,R4Au 0). 

Then we deduce that for any \y, h) e Bidzo) x [Ar/R , (A,, + 1)R ), there exists a time f G 
(-yS2/?^,0) such that 

h = inf Wy = Wy [X, 1) , 

B„,Kfco)x(-Q-2R^'] ' ■ V / 

where the infimum is achieved at an interior point x e Bp^Rizo). This means that (x't) is a 
parabolic contact point, i.e., (x, G yi. According to Lemma 4.2 we observe that 

y = exp^ -7?^ Vm (x, f) , and xiCvA(y). 

Now, we define the map (p : ^q.,r,„,r2(zo, 0) — > M as 

^{x, f) :- exp^ [ -R^Vu{x, t)\, 

and the map O : ^a,R,a,2R2(zo, 0) ^ M x R as 

OCx, f) := [(;*(x, t),--cf (x, 0(x, f)) - TzR^uix, f) + C,,f 



We also define 

:A := I (x, t) G Ki3^R,p,R2(zo, 0) : 3ye 5r(zo) s.t. w/x, f) = inf Wy < (A;, + 1)^M . 

(^ B„,«(zo)x(-ff2«-,'] ' j 

According to the argument above, we have proved that for any (y, s) e Br(zo)x(-(A,j + 1)R^, -A,jR^), 
there exists a point (x, f) G ^ such that {y, s) - <l>(x, t), that is, 

Br(zo) X (-(A, + 1)R\ -A,R^) c (D(^). 

Thus, the area formula provides 

(24) \Br(zq)\ R^ < f -H" [^ n d)"' (y, i)l t/VCy, i) = r Jac <D(x, O^V(x, 0- 
We note that 

(25) ^ c y[ n /:^,R,;j2«^fa), 0) n {m < m,) 

for M,, := 2(A;, + 1) since \R^u(x, t) < Wy(x, t) < (A,, + \)R^ for (x, f) G :R. 
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Next, we claim that for (x, t) e y[ 
(26) 

1 



From Lemma [43J if (x, f) e J[, then we have 

^2 _ 

— Vm(x, = -(iy(x)Vrfv(.x) for _v := 4>{x, t) e Br{zo); x i Cut(y), 



n+I 



and hence 



(27) 



-|VM(x,f)| - dy(x) < d(y,zo) + d(zo,x) < R + aiR < laiRij. 



Using Lemma 4.4 (with A = A/n) and pT] ), we deduce that for (x, f) e y[, 
(n + l)Jac<I)(x, 0"+" 



<S^ 



(n- 1)kR^\Vu\'' 



(riR^ (A ^ \ 



6m ,^f l(n-l)KR^\Vu\] 



nR^ I A 



6« 



< ^ {lai V^/?()) i — [- AM - 5,m] + -^ + n,yf (lai ^/kRqJ 



since ^(t) and ^(t) are nondecreasing for t > 0. This proves ( [26) . 
Lastly, we shall show that for (x, t) € J[, 



(28) 



-AM < M"(D2m) + -V^ (2a'i -y^/^o) ■ 



2.7 



Indeed, for (x, f) e J[, we recall Lemma 4.2 again to see 

D^[u+ —dl\{x,t)>0 for3;:=0(x,f); x ^ Cut(3;), 

i.e., the Hessian of R^u + d^. at (x, t) is positive semidefinite. From Lemma 
follows that 

D\{x,t) > -^D^i}-dl\{x) > -^j^[yrKdy{x)) g, > -^Jf{2ai ^/^Ro)g, 

Let fj.1 be the largest eigenvalue of D^u{x, t). If jUi > 0, then we have 

2 
NV{D^u{x,t)) > Am - (n - l)K—M'{2ai yf^^Ro) 



and d27l, it 



R2 

A 2 

> -Au-riA—rJ^ilaiy/KRo). 
n R^ ^ ' 



If yUi < 0, then we have 



2 1 o 

Mr{D^u(x,tf) = Aam > -nK—.ye{lax ^kR(^ > - l\u - nK—.^e {lax V^^o) , 



which proves (|28]l for (x, t) e y{. Therefore, the ABP-Krylov-Tso type estimate ( [23| l fol- 
lows from (|24ll, dBll (|26ll and (|28ll. n 
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We modify the barrier function of [WJ (see IKKLI ') to construct a barrier function in the 
Riemannian setting. First, we fix some constants that will be used frequently (see Figure 
n\; for a given < 77 < 1, 

11 n"* 9 

ai :- — , 0-2 := 4 + 77^ H , /3i := - and /32 := 4 + 77^. 

77 4 T] 

Lemma 4.6. Assume that 

Sec > -K on M, for k>0. 

Let Ro > and < 77 < 1. For zo e M, and < R < Rq, there exists a continuous function 
v,^{x,t) in K^^K a^giizo,/^!^^), which is smooth in (M\Cut(zo)) n K„^i(„^j(i{zQ,li2R^), such 
that 

(a) v^(x, f) > in /:„,«_„,R2(zo,y62/?^) \ Kp.R^pMzo^PiRh 

(b) v^(x,t)<OinK2R(zo,/32Rh 

(c) R^ [M^(D\) - 5,v,}+p+2(77+l)A^(2ai ^/J^Rq) <Oa.e. m /:^,«,^,«2(zo,y62-R')\^2R(z„, '^R^), 

(d) R^ [M^(D\) - d,v^] < C, a.e. in Kp^R^pMzo^PiRh 

(e) v^(x,t) > -C^ in ^„,R,ff,fi2(zo,j02^^), 

where J^(t) — t coth(T), and the constant C,, > depends only on rj, n. A, A, ^fKRo (inde- 
pendent ofR and zo)- 

Proof As in iTWl Lemma 3.22] and IKKLI Lemma 4.1], we consider 

h(s, t) := -Ae-"" 1-4 .. \ ;, exp (-a-) for t > 0, 

and define 

^,,(.,0 - h(s,t) + Ct in [0,yS2] X [0,/?2]\[0, 'i] X [0, ^], 

where C :- 12/77^ + 2(n + 1)AJ^(2q'i V^/?oj, and the positive constants A,m,l,a (de- 
pending only on 77, n. A, A, V'i'^o) will be chosen later. Extending ip,j smoothly in [0, a^] x 

4 
[-x.jSi] to satisfy (a) and (e), we define 



4 

dl(x) t ^ 



v^ixj) :=<^,, 



for (x, t) e /:„,«_ „^R2(zo,y82^), 



R^ 'R^j 

where d,„ is the distance function to zq- We may assume that (fiijis, t) is nondecreasing with 

respect to s in [0, aj] x [-'■j,/32]- 
We recall that 

(ZJ2 [dl/2) (x) ■ ^, ^) = {d,„D%^(x) ■^,^) + (Vflf,„(x), ^f , V^ € r,M, X i Cut(zo), 

and 

M^ (£»' « /2) (x)) < tjA^ ( ^/^d,, (x)) < nKM' (ai V^T^o) , Vx e B^,r(zo)\ Cut(zo), 



from Lemma[Z7| Following the proof of BKKLI Lemma 4.1], and using Lemma[2.15 (a). 



we can select positive constants A,m,l,a, depending only on ri,n,A,A, -s/kRo, such that 
(b), (c), and (d) hold. For the details, we refer to the proof of flKKLI Lemma 4.1] (see also 
llWl Lemma 3.22]). n 

We approximate the barrier function v,, by a sequence of smooth functions as Cabre's 
approach in BCal since v^(x, t) - ip^l ^Ij— > ^ I is not smooth on Cut(z()). We note that the 
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cut locus of zo is closed and has measure zero. It is not hard to verify the following lemma 
and we refer to BCal Lemmas 5.3, 5.4]. 

Lemma 4.7. Let zo e M, R > and let (p : R^ x[0,T] -^ R be a smooth function such 
that (p(s, t) is nondecreasing with respect to s for any t € [0, T\. Let v{x, t) :— (f{(£ (x), n. 
Then there exist a smooth function ijj on M satisfying 

< i/r < 1 on M, ifr = I in Ba k{zo), and suppi/r C BwJzo), 

and a sequence [wk]'^^^ of smooth functions in M X [0, T] such that 

Wk — » i/'v uniformly in M X [0, T], 

dfWk — » i//dtV uniformly in M X [0, T], 

D^Wk < Cg, in M X [0, T], 

D^wt -^ D\ a.e. in B/j,fi(zo) X [0, T], 

where the constant C > is independent ofk. 



The following proposition is obtained by applying Lemma 4.5 to m + v^, with v^, con- 
structed in Lemma 4.6 and translated in time, due to the approximation lemma above. 

Proposition 4.8. Assume that 

Sec > -K on M, for k > 0, 

and that F satisfies (HI). Let < rj < \ and ^q.,r, 0.2*2(20, 4/?^) c KrJxq, to) <z M x R. Let 
u be a smooth function on K^^k,^^_k2{zo,4-R^) such that 

F{D\)-d,u<f in Kp^R,p,R2{zoAR\ 



and 



u>Q in K^,R,a.MzoAR^)\Kp,R,f},R2(xoAR^), 



inf u < I. 

Kir(z«AR2) 



Then, there exist uniform constants M,^ > 1, < yU^, < 1, and < e^, < 1 such that 

|{M<M,}n/:,R(zo,0)| 



KaR,AR2(zo,4/?2)| 



■^i^l^ 



provided 



i 



\fi\R^r 



infl+l 



<e„; r:=max(/,0). 



where 6 :- \ + log2 cosh(4 ^/kRo), and M^ > 0, < yu,,, e,j < I depend only on rj, n. A, A 
and yfi<Ro. 



Proof Let v^ be the barrier function as in Lemma 4.6 after translation in time (by -rj R ) 

We 



4.7 



and let {w,i)^j be a sequence of smooth functions approximating v,, from Lemma 

notice that m + v^, > in Ka^R^a^R2(zoAR^)\Kfj,R,fi^R2(zt_)AR^) and inf (m + V;,) < L We 



K2r(zoAR'-) 



4.3], and use the dominated convergence theorem to let k go to +00 due to Lemma 4.7 
Thus we obtain 



can apply Lemma 4.5 to u + wt after a slight modification as in the proof of MKKLI Lemma 
dom: 

/ 

Jtu-l-v 



\Br(zo)\ ■R^<Ci 



-i'„<M„)nif^j„^^„2(zo,4«-) 



[r^ [MriD^u + D^V;,) - d,{u + V;,)} + Ci] 
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where Ci := J5^"+' (2 V^/?o) l{2A)"^\ and Ci := lllrf + 2{n + \)KJ^ {l yfl<RQ) . Using 
Lemma 2.15 (HI), and the properties (c),(£/) of V;, in Lemma 4.6 we have 

C V "itt+l 

|Br(zo)I ■R'<C, [r^ [M-{D\) - d,u} + R^ [MHd\,) - drv,) + Cz}^ 



Ci r \{r^ [piD^u) - d,u} + R^ {m^(D\) - d,v^ + Czl 

"I 



\R^r + (C, + C2)XE,r'^ 



where Ey := {m + v^, < M,,) n (^/j,R,/j,fi2(zo,47;2)\7i:^s(2o,0)) and £2 := {m + V;, < M,,) n 
K,jr{zo, 0). Then, it follows that 

\BRizo)\ ■ R^ 






i 



< C3 1 t IfiR'f 



i2D2^+|«9+l 



+ c, 



1^2 



K/S,R,/3,R2(Z0,4/?2)|" 



for :- I + log2 cosh(4 V^/^o) > 1, where a uniform constant C3 > depending only on 
T], n. A, A and ^]KRQ may change from line to line. Therefore, Bishop-Gromov's Theorem 
|2.3|implies that 



\E2 



Kp^R,pMzoAR^) 



i 



'p,R.P2R2{Z,^R^) 



\fi\R^r 



nfi+I 



1 |gfi(z())l ■ R^ 
Ci\Kp^R,^_R2{zaARH 



11/1 \'°°^^ 1 



Pi 



CiDXP 



--: 2<^' 



for £) := 2" cosh" '(2 ^fKR^)). By selecting e,, := /v^^'"*' , we conclude that 



^''l 



\[u + V, < M,} n K,jR(zi.u 0)| |{m < M^} n K,^r{zo, 0)| 

|X'/j,R,/3,R2(zo,4/?2)| 



^B,«,ft«Kzo,47;2) 



for M^ :- Mri+Cri depending only on 77, n. A, A and ^]KRQ since V;, > -C^ in Kp^R p^RiizQ, A-R^) 
from Lemma l46l 



4.2. Parabolic Harnack inequalities for viscosity solutions. Now we shall prove Propo- 
1 4.9 which is a key ingredient in proving Theorems |1.1| and 1.2 



sition 



Proposition 4.9. Assume that 

Sec > -K on M, for k>0, 

and that F satisfies (HI). Let < 77 < 1 and K^^r g,^_R2{zQ,4R^) C KR^ixQ, to) C M X R. For 
f & C (A's„(xo, ?())) , let u e C {Kr^(xo, to)) be a viscosity supersolution of 

F{D\)-d,u^f in K,^r^,^r2(zoAR\ 



such that 



u>0 in K,,R,a2R<ZoAR^)\Kp,R,p^R2(zoAR\ 
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and 



inf M < 1. 



Then, tliere exist uniform constants M^ > 1, < yu^, < 1, and < e^, < 1 such that 

|{M<M,}n/:,fi(zo,0)| 



provided 



(29) 1 1 \Rirf'' 



[i 



KRAXQ^tQ) 



■f^n' 



^^v. 



where 9 :- \ + logj cosh(4 yfiiRo), and M^ > 0, < //,,, e,j < 1 depend only on rj, n. A, A 
and yfKRo. 

Proof. It suffices to prove the proposition for F - Mr from (HI) (or (HI')). Setting 

a\ :- (ai + /3])/2, and q'2 :- (ca +j62)/2, we define 

Q X (ro, Tz] := K^,r,^^r2(zoAR\ and // x (r,, Tj] := K^^R^/sMzo^'^^^y 
We note that u and / belong to C (Q x [To, Tz]) , and we denote by oj the modulus of 

continuity of m on Q x [Tq, T2], which is nondecreasing with (jj(0+) = 0. 

For e > 0, let u^ be the inf-convolution of u with respect to Q x (To, T2] as in (j9|. 
According to Lemma 3.6 there exists eo > such that if < e < eo, then u^ satisfies 

M'(D^Us) - d,u^ < fe a.e. in K/^^R/j^niizoAR^), 

where /^ is defined as follows: for (x, f) e Kp^j^ p^gi{za,AR^), 

fs(x,t):^_ sup f + 2nAK(xJ{2^/msy, m:=\\u\\,- ,, ^j^i^V 

and we recall that A1" is intrinsically uniformly continuous with respect to x with a>M- = 0. 
Using ([8| and (|29|, we have that 

<2(|[-',,=:l„ 

and hence for small e > 0, 
(30) 

< 2e,„ 

since fs converges uniformly to / in Kp^j^ p^jfiizaAR^)- For a fixed (5 > 0, we may assume 
that for small e > 0, 

u,>-5 in Ki,^R^^^K^{zoAR^)\Kp,R^ii2R^{zi)AR^), 

and 

inf Ug, <\ + S 

KiRizoAR-) 

■?2\ 



since Mg converges uniformly to u in K^^f^^^^KiizoAR ) from Lemma 



3.1 
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Now, we fix a small e > 0. According to Lemma 3.2 (c), there is a smooth function ip 
on M X (-00, Ti] satisfying < (A < 1 on M x (-oo, Ta], 

i/r=l mKp^Rfi^gi(z()AR^), and suppi/r c /:3,R,a^jR2(zo,4/?^), 

and we find a sequence {w^tl^i ^^ smooth functions on M x (-oo, Ta] satisfying 



IVwtl + |5,wt| < C 
D^wt < Cg 



uniformly in M x (-oo, T2\ as A: — > +oo, 

inMx(-co,r2], 

a.e. in Kp^R^p^j(i{zoAR^) as ^ -> +oo, 

inMx(-oo,r2], 

a.e. in Kp^gp^jiiizoAR^) as fc -> +oo. 



where the constant C > is independent of k. For large k, we may assume that 



Wk > -26 mKa^R^aMzoAR')\Kp,R,p,R2{zoAR'), 



inf 

KMzoAR^) 1 + 45 



< 1, 



and 



i 



<4e;„ 



to 



1 +45 



(for large k) to 



le^/;^ {m-(z)2w,) - d,wk 

where we used the dominated convergence theorem to obt ain t he last estimate from pO| 

Wk + 26 
Selecting 6;, > small enough, we apply Proposition 4.P '" 

obtain 

|{w, + 25<(l+45)M,}n/i:„R(zo,0)| 

\Kp,R,/sM^oAR^ 
By letting k — > +oo, we have 

\[u, + 6 < {I + 46)M^} n K^r(zo,0)\ 

\K/s,R,fl,R2{zoAR^ 

Since m^ converges uniformly to u in Kq^r^,r2(zoAR^), we let s 
Bishop-Gromov's Theorem 2.3 to deduce that 

\[u<M^]nK,jR{zo,0) 



■l^i- 



^i,r 



and 5 — > 0, and use 



\Ka,R,a2R-iZQ'^R') 






for D :-2" cosh" (4 -s/kRo), which finishes the proof. n 

Sketch of proof of Theorems |l.l| and [L2| 

Theorems |1.1| and |1.2| follow from Proposition 4.9 and a standard covering argument 
using Bishop and Gromov's Theorem 2.3 Indeed, we first prove a decay estimate for the 
distribution function of a viscosity supersolution u to F{D^u) - d,u = / in K2r{x(iAR^) by 
using Proposition 4.9 and a parabolic version of the Calderon-Zygmund decomposition in 
BChL Then, the weak Harnack inequality in Theorem [L2] easily follows from Bishop and 
Gromov's Theorem 2.3 To complete the proof of Theorem |1.1[ we apply Proposition 4.9 
and obtain the same decay estimate for w :- C\- C2U (for Ci, C2 > 0), which satisfies 

M'(D^w) - drw = -C2 {M^(D^u) - d,u] < -Ci {f(D^u) - d,u] = -Cif 

in the viscosity sense. For the detailed proofs, we refer to flKKLIIWll . 
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5. Elliptic Harnack inequality 

Using the sup and inf- convolutions, we shall prove Harnack inequality of continuous 
viscosity solutions to elliptic equations from a priori estimates. We recall viscosity solu- 
tions for uniformly elliptic operators. 

Definition 5.1 (Viscosity sub and super-differentials, fAFSl). Let Q. a M be open and let 
u : Q. ^t Wbe a lower semi-continuous function. We define the second order subjet ofu at 
X e Clby 

J^''u(x) := {(v<^(x),dVw) 6 r,M X SymTM^- : ip € C^(Q), 
u — ip has a local minimum at x) . 

If ((,A) e Sr^'~u(x), ^ and A are called a first order subdifferential and a second order 
subdijferential ofu at x, respectively. 

Similarly, for a upper semi-continuous function m : Q — > R, we define the second order 
superjet of u at x & Q.by 

J^-^u(x) := ((Vy)(jc),DVw) e T_,M x SymTM., : (f e C^(Q), 
u — <p has a local maximum at x] . 

We quote the following local characterization of J'^'^u from BAFSI Proposition 2.2]. 

Lemma 5.2. Let u : Q. ^ 'Rbe a lower semi-continuous function and x e M. The following 
statements are equivalent: 
(a) (^,A) e J^--u(x). 

Definition 5.3 (Viscosity solution). LetF : MxRx TM x Sym TM -^ R and let D. c M 
be an open set. We say that a upper semi-continuous function m : Q — > R /i a viscosity 
subsolution of the equation F(x, u, Vm, D^u) — on Q. if 

F(x,u(x),(,A)>0 

for any x e Q. and (^, A) e J'^'^u(x). Similarly, a lower semi-continuous function u : Q. —> 
R is said to be a viscosity supersolution of the equation F{x, u, Vm, D u) = Q on Q. if 

F(x,u(x),(,A) <0 

for any ^ e Q and (^, A) e J'^'^u(x). We say u is a viscosity solution ifu is both a viscosity 
subsolution and a viscosity supersolution. 

As in the parabolic case, we use the sup and inf- convolution to approximate continuous 
viscosity solutions. Let Q c M be a bounded open set, and m be a continuous function 
on Q. For e > 0, we define the inf -convolution of u (with respect to Q), denoted by u^, as 
follows: for xq € Q., 

Ui-ixo) := inf \ u(y) -i- —d^(y,xo) \ 
If yo is a point to realize the above infimum, then we have 

1 2 1 2 

Ue(xq) = uiyo) -H —d (y(,,XQ) < u(y) H- —d (y,XQ) Vy e D., 
2s Is 

which means u has a touching paraboloid -^(i^(y, xq) -h uiyo) -i- j^d^(yQ,X()) at yo from 
below. 

Now we state the elliptic analogue of the results in Section[3]without proof. 



(b) u (exp,^) > u(x) + <^,^) + i(A^,^) + o(|^|2) as T,M 3 ^ -^ 0. 
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Lemma 5.4. For u e C (Oj , let u^ be the inf-convolution of u with respect to Q and let 
X(j e Q. 

(a) IfO < E < s' , then Me'(ji;o) < MeCxq) < u{xq). 

(b) There exists yo € Q. such that u^ixQ) - uijo) + i^d^ijo, xq). 

(c) d'^(yQ,xo) < 2s\u(xo)- u(yo)\ < 4e||M||L».(n). 

(d) Ug 1 u uniformly in Q. 

(e) Ue is Lipschitz continuous in Q.: for xq, x\ e Q, 

3 

\uc(xq) - Us{x\)\ < —d\dirt\{D)d{xQ,x\). 

2e 

Lemma 5.5. Assume that 

Sec > -K on M, for k > 0. 

For u e C (O.) , let u^ be the inf convolution of u with respect to Q, where Q. d M is a 
bounded open set. 

(a) Ug is semi-concave in Q. Moreover, for almost every x € H, u^ is dijferentiable at x, 
and there exists the Hessian D^Ue{x) (in the sense of Aleksandrov-Bangert's Theorem 
\2.8\ such that 

u, (exp, ^) = u,{x) + <Vm,(x), ^) + ^-{A{x) ■^,^)+o (l^l^) 
as^ e T^M -^ 0. 

(b) D Ue(x) < - V^diam(Q)coth( V'fdiam(Q)Wj a.e. in D.. 

(c) Let H be an open set such that // C Q. Then, there exist a smooth function t// on M 
satisfying 

< t// < I on M, t// = I in H and supp i// C Q., 
and a sequence {wk}'^^^ of smooth functions on M such that 

Wk — > ip^E uniformly in M as k -^ +co, 

\Vwk\ < C in M, 

D^Wk < Cg in M, 

D^Wk —> D^Ug a.e. in H as k —> +oo, 

where the constant C > is independent ofk. 

Proposition 5.6. Assume that 

Sec > -K on M, for k > 0. 

Let H <Z M be a bounded open set such that H d Q.. Let u e C ( Q j , and let oj denote 

a modulus of continuity of u on Q, which is nondecreasing on (0, +oo) with w(0+) = 0. 
For e > 0, let u^ be the inf-convolution of u with respect to Q. Then, there exists sq > 
depending only on \\u\\^„/^\,fi^ andQ., such that if Q < e < sq, then the following statements 

hold: Let xq e H, and let yo & Q. satisfy 

1 2 

u^ixo) = M(yo) + —d {yQ,XQ). 
2e 

(a) We have that yo € Q, and there is a unique minimizing geodesic joining xq to ya. 

(b) If{(,A) e JT^'^MeCjcq), then we have 

yo = exp (-eO- 
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(c) If{(,A) € ^ Usi^o), then we have 

where L^cj-o stands for the parallel transport along the unique minimizing geodesic 
joining xq to yo = exp (-ef ). 



Lemma 5.7. Under the same assumption as Proposition 5.6 we also assume that F satis- 
fies (HI) and ( |H2| l. For f & C (Q), let u e C (O.) be a viscosity supersolution of 

F(D^u) = / in Q. 
IfO < E < £(), then the inf convolution Ug (with respect to Q.) is a viscosity supersolution of 

FiD^Ui,) = f,, on H, 
where sq > is the constant as in Proposition \5. 6\ and 



feix):— sup / + Dp (2 y/ms) + 2nAK co (2 y/ms) ; m := ||m||^„/^\. 
Moreover, we have 






F(D Ue) < fs a.e. in H. 



In particular, Lemma 5.7 holds for Pucci's extremal operators according to Remark 3.5 



The following proposition is quoted from the proof of I WZ[ Proposition 4. 1], which is 
a main ingredient in the proof of a priori Harnack estimate. 

Proposition 5.8. Assume that 

Sec > -K on M, for /c > 0, 

and F satisfies (HI). Let u be a smooth function satisfying FiD^u) < f in B2r{xq) C M. 
Then, there exist uniform constants Mq > 0,0 < /io < 1, and < eo < 1, depending only 
on n, A, A and yflcR, such that if for any Byizo) C B2r(xq), 



and 



then we have 



M > in BJzo), inf m < 1, 

B,/2(zo) 



|rVl"1 <eo; r:=max(/,0). 



B3rll(Zo) 

|{M<Mo)nB,/g(zo) 



:yUo, 



\BriZ0)\ 

where 9 :- \ + logj cosh (8 yf^Pi] . 

Lemma 5.9. Assume that 

Sec > -K on M, for /c > 0, 

and F satisfies (HI). For f & C (B2r(xo)) , let u e C (B2r(xq)) be a viscosity supersolution 
of 

F(D^u) = / in B2r(xo). 

Then, there exist uniform constants Mq > 0,0 < ;Uo < 1^ ond < eo < 1; depending only 
on n. A, A and ^/kR, such that if for any Brizo) C B2r(xq), 

u > in B,.(zq), and inf u < I, 

Briiiza) 



then we have 



provided that 
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\{u < Mo) n B,/8(Zo) 



\Br(ZQ)\ 



■-1^0, 



i 



B2r{xq) 



<eo; r:=max(/,0), 

where 6 :- \ + logj cosh (8 ^[i<.R\ . 

Proof. It suffices to prove the proposition for F = A1" from (HI). Set 

Q. := Bvr/sfa)), and H := £3^/4(20)- 
We note that u and / belong to C (Q j , and denote by w the modulus of continuity of u on 

Q, which is nondecreasing with w(0+) = 0. 

For small s > 0, let u^ be the inf-convolution of u with respect to Q. According to 
Lemma [5^ we find eq > such that for < e < so, m^ satisfies 

M'iD^Ue) < fe a.e. in £3^/4(20), 

where /e is defined as follows: for x e £3^/4(20), 

fi.(x) := _sup / + 2nKK(x)(2 yfms^ ; m := I|m|Iz.»(b^^^^(,^)), 



BiJi^M 



and we recall that Al is intrinsically uniformly continuous with respect to x with ojm- = 0. 
Using (|7|, we have that 

2 — — 

(t) ff |.rr)"<2-4/?2(/ |/+|«e]'% 860; 0:=l+log2Cosh(8V^7?), 

\ ^ I \Jb_v/4(Z0) / \JB2r(xo) I 

and hence for small e > 0, 

since fg converges uniformly to / in B3,/4(z()). For a fixed 5 > 0, we may assume that for 
small e > 0, 

Ue ^ -S in Bjri&izo), and inf Ue < I + 6 

Br/l(Zo) 



since u^ converges uniformly to u in Bjr/sizu) from Lemma 5.4 



Now, we fix a small e > 0. According to Lemma [53] (c), there is a smooth function tf/ 
such that < iff < I on M, 

iff si in B3,./4(zo), and suppt// c Bj^/aizo), 

and we approximate tfrug; by smooth functions Wk on M satisfying 

Wk -^ 4'^e uniformly in M as A: — » +00, 

\Vwk\ <C in M, 

D^Wk < Cg in M, 

D^Wk -^ D^Ue a.e. in B3,/4(zo) as A: — > +00, 

where the constant C > is independent of k. For large k, we may assume that 
Wk > -2(5 in Brizo) and inf Wk < I + 26, 

Br/2(Z0) 
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and 

\2 



f?)(f |{M-(^v)}r)% 10.0, 



where we used the dominated convergence theorem to ob tain the last estimate from ( [3T| l. 

Wk + 26 

By selecting eo > small enough, we apply Theorem 5.8 to (for large k) to obtain 

1+45 



that 

|{w^ + 2(?<(l+45)Mo)nfi,/8(zo)| 

|B.(zo)l " ^°' 

By letting ^ — » +oo, we have 

|{M, + 5<(l+45)Mo)nB,/8(zo)| 

Since m^ converges uniformly to u in BT,-/fi,{zo), we let e ^ and 5 — > to conclude that 

\{u < Mq) n g,/8(Zo)| 

\Br{Z0)\ -^° 

which finishes the proof n 

Proof of Theorems 11.31 and 11.41 Theorems 11.31 and 11.41 follow from Lemma 15.91 and 
a standard covering argument using Bishop and Gromov's Theorem |2. 3 1 (see ICal for in- 
stance). D 
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